Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 October 2008 (MN KTeX style file v2.2) 



Stellar disruption by a supermassive black hole: is the light curve 
really proportional to r^^^l 



G. Lodato\ A. R. King^ and J. E. Pringle^'^ 

' Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH 
^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OHA 



ABSTRACT 

In this paper we revisit the arguments for the basis of the time evolution of the flares ex- 
pected to arise when a star is disrupted by a supermassive black hole. We present a simple 
analytic model relating the lightcurve to the internal density structure of the star. We thus 
show that the standard lightcurve proportional to r^^^ only holds at late times. Close to the 
peak luminosity the lightcurve is shallower, deviating more strongly from r^^^ for more cen- 
trally concentrated (e.g. solar-type) stars. We test our model numerically by simulating the 
tidal disruption of several stellar models, described by simple polytropic spheres with index 
y. The simulations agree with the analytical model given two considerations. First, the stars 
are somewhat inflated on reaching pericentre because of the effective reduction of gravity in 
the tidal field of the black hole. This is well described by a homologous expansion by a factor 
which becomes smaller as the polytropic index becomes larger. Second, for large polytropic 
indices wings appear in the tails of the energy distribution, indicating that some material is 
pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, 
the r^^^ lightcurve is achieved only at late stages. In particular we predict that for solar type 
stars, this happens only after the luminosity has dropped by at least two magnitudes from the 
peak. We discuss our results in the light of recent observations of flares in otherwise quiescent 
galaxies and note the dependence of these results on further parameters, such as the star/hole 
mass ratio and the stellar orbit 
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1 INTRODUCTION 

X-ray flares from quiescent (non-AGN) galaxies are often inter- 
preted as arising from the tidal disruption of stars as they get close 
to a dormant supermassive black hole (SMBH) in the centre of the 
> galaxy ( Komossa & Bade 1999). Similar processes also occur on 
' much smaller scales, such as in c ompact binary systems , where the 
black hole is only of stellar mass jRosswog et aklboOSiffl) . 

The pionee ring work by iLacv et al.l h982h : iReesI ( Il988h : 
PhinnevI 1 19891) a nd, fr om a numerical point of view, by 
Evans & Kochanekl ( [198^ have set the theoretical standard for the 
interpretation of such events. In particular, a distinctive feature of 
this theory is an apparent prediction for the time dependence of the 
light curve of such events, in the form Lit) oc r''^ (note that the 
original paper by|Rees 1988 quotes a /"''^ dependence, later cor- 
rected to r^'^ by Phinnev 1989). Since then, a r^'^ light curve is 
generally fitted to the observed luminosities of events interpreted 
as stellar disruptions. In this paper, we revisit the theoretical ar- 
guments behind_such scaling and we show (as also originally ar- 
gued bv lReesiri988h that the light curve does not need to have this 
scaling and, in particular, that it critically depends on the internal 
structure of the star being disrupted. We provide a simple model 
to calculate the light curve starting from the density profile of the 



star and we show that more centrally concentrated stars tend to pro- 
duce shallower light curves. We further supplement our model by a 
numerical calculation of the process, using Smoothed Particle Hy- 
drodynamics (SPH). 

W e start by briefly summarizing the main argument of iReej 
( Il988b . Let us consider a star, originally in hydrostatic equilibrium 
at a large distance from the black hole. Since pressure forces and 
the internal self-gravity of the star are in equilibrium, the only un- 
balanced force is the gravitational pull of the black hole. The vari- 
ous fluid elements of the star therefore move in essentially Keple- 
rian orbits around the black hole, each one with its own eccentricity, 
that is initially very close to the eccentricity of the centre of mass 
of the star (in the following, for simplicity, we make the simple as- 
sumption that the centre of mass of the star is in a parabolic orbit 
around the black hole). Therefore, the distribution of specific me- 
chanical energy within the star is very narrow around the energy of 
the centre of mass. As the star moves closer to the black hole, the 
various Keplerian orbits tend to be squeezed, perturbing the hydro- 
static balance. Pressure forces then redistribute energy inside the 
star, therefore widening the specific energy distribution. After the 
encounter, the star is thus characterized by a much wider distribu- 
tion of internal energies, with part of the fluid having a negative en- 
ergy (and therefore being bound to the black hole) and part having 
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a positive energy (and therefore remaining unbound). In the picture 
of|Rees ( 1988) it is this energy distribution (and only this) that de- 
termines the light curve of the event. Indeed, after the encounter the 
fluid elements again move in Keplerian orbits (but with their new 
energy). The bound elements then come back close to pericentre 
after a Keplerian period T, linked to their (negative) energy E by: 



^ l/ 2;rGMh \ 



2/3 



(1) 



where Mj, is the black hole mass. The mass distribution with spe- 
cific energy dM/d£ then translates, through Eq. ([T]l, into a mass 
distribution of return times dM/dr. The next fundamental assump- 
tion is that once the bound material has come back to the pericen- 
tre it loses its energy and angular momentum on a timescale much 
shorter than T, thus suddenly accreting onto the SMBH and giv- 
ing rise to the flare. The mass distribution of return times is there- 
fore effectively the mass accretion rate of the black hole during the 
event, from which the luminosity can be easily computed. We thus 
have: 




ToBH 



AM _ AM dE 
dT ~ dE dT 



OnGM^dM 
3 dE 



(2) 



In order to obtain the 'standard' r^^^ light curve, we then have to 
make the second fundam ental assump tion that the energy distribu- 
tion is uniform. Note thatlReesl( ll988l) did not show that this should 
be the case, and only assume d it for s i mpl icity. Later, the numeri- 
cal simulations by Evans & K ochane^ |l989) apparently showed a 
uniform energy distribution, hence suggesting that the light curve 
is generally proportional to In the following, we first show 
analytically that the energy distribution need not be uniform, but 
depends on the properties of the star, and in particular on its inter- 
nal structure. We then show numerically that in fact it is not uni- 
form and does depend on the properties of the stars, in a way that 

approximately reproduces the analytical results. 

St arting from the pioneering work of ICarter & Luminell 

( Il982l. [ 19830 . numerical simulations of this process have been 



performed in a variety of studies jBickne 1 & Gingold 
Evans & Kochanek'l lgSa : iLaguna et al.lll993al lbl: lAval et al.l 



1983 



200C 



Bogdanovic et al. 2004D. Mo st of these studies used Smoothed Par- 
ticle Hydrodynamics (SPH. 'Mona gharj 1 1992h . mostly because of 
its ability of following the system over a wide range of physical 
scales, where much of the simulated region is essentially 'empty'. 
These various attempts have considered t he effects of vary- 
ing the orbital parameters of th e encounter jBicknell & Gingol3 
1 19831 : 1 Evans & Kochanekl [T989h , of the inclusion of relativistic 
terms in the equation of motion iLaguna et alJI 19933 15 1 Aval et al.l 
I2OOOI) and have describ ed the expected observational outcome 
I Bogdanovic et al.l l2004h . However, surprisingly, no attempt has 
been made at exploring the effect of varying the internal structure of 
the star. Indeed, all such analyses have considered a simple poly- 
tropic model for the star, with 7 invariably set to 5/3 (although 
iRosswog et al.l2008al lblhas consider the encounter of a white dwarf 
with a stellar mass black hole, which is in the much smaller mass 
ratio regime). 

The paper is organized as follows. In section 2 we describe 
our analytical model to derive the energy distribution of the dis- 
rupted debris. In section 3 we describe our numerical code and the 
set up of our simulations. In section 4 we describe the results of 
the simulations. In section 5 we discuss our results and draw our 
conclusions. 



Figure 1. Schematic view of the geometry of the system. The radius of the 
star is Rj,Jhe SMBH is on the right, at a distance Sp 2> S*. 



2 THE PROCESS OF TIDAL DISRUPTION OF A STAR 
BY A SMBH 

A simple and instructive way to consider the process is by treating 
the interaction of the star with the black hole under the impulse ap- 
proximation, that is assuming that the interaction occurs in a very 
small time span as the star gets close to pericentre. This approx- 
imation is probably appropriate for highly hyperbolic encounters, 
but is only approximate for the parabolic case considered here. We 
therefore expect the details to differ somewhat from what is de- 
rived here below, even if the qualitative behaviour is correct. In 
this approximation, the motion of the star is simply a straight line 
until it reaches pericentre, at which point it is subject to a short 
impulse that deflects the various fluid elements, each of them in- 
dividually conserving their specific energy. Until it reaches peri- 
centre, therefore, the structure of the star is essentially unchanged 
and it keeps its original radial density profile and its initial radius. 
The key thing to realize here is that the spread of specific energy 
the star has reached just before the impulse is simply given by the 
different depths at which the various fluid elements are within the 
black hole potential well, because they all share the same veloc- 
ity and therefore the same kinetic energy. This was c learly realized 
by iLacv et aP ( Il982h and lEvans & KochaneS l ll989h who pointed 
out that "The spread in specific energy of the gas... is given by 
the change in the black hole potential across a stellar radius". The 
iinpulse occurs istantaneously and therefore does not modify the 
kinetic energy of the fluid emlements but imparts some degree of 
rotation in the star , with » (2GMh/j?p)'^^ , where Rp is the peri- 
center distance (cf. lEvans & Kochanekl ^ 1989*). We consider here the 
case in which the stellar radius R^, is much smaller than Rp, which 
is appropriate for the case of stellar disruption by a supermassive 
black holes (although not for the case of compact binaries). This 
allows us to easily estimate the expected energy spread, as 



AE ■ 



d£p 
dr 



(3) 
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Figure 2. Left: Distribution of internal energy for polytropic stars with different indices. Solid line: y = 5/3, long-dashed line: y = 1.4, short-dashed line: 
y = 4/3. Center: corresponding evolution of the accretion rate for the same three cases. The red line indicates for comparison a simple f"'''^ power-law. Right: 
time evolution of the power-law index n = dlnm/dln t. As can be seen the value n = -5/3 is only approached at late times. 



where iSp = GM^^Ir is the potential energy due to the black hole and 
Ar max = R* is the ma ximum deviation from the pericentre distance 
(cf. lLacvetalJl982l) . We thus expect the energy distribution to ex- 
tend roughly between -A£ and Aii. This simple result has been 
obtained in all the early analyses of the problem. However, in the 
same approximation as before, we are also able to derive the whole 
energy distribution starting from the density, by calculating what 
is the fraction of stellar mass at a given Ar from the centre. Fig.[T] 
illustrates the geometry. It can be easily shown that 



— = In p{r)rdr, 
dAr 



(4) 



where p(r) is the spherically symmetric mass density of the star. 
The relation between the distribution of Ar and the distribution of 
energy E is simply given by: 



dM 
d£ 



dM R^, 
dA^AE' 



(5) 



where AE is given by Eq. (O. 

It is useful to introduce dimensionless quantities. We then de- 
fine 6 = -E/AE (where we have also included a minus sign be- 
cause we are interested in material with negative specific energy) 
as our dimensionless energy, x = Ar/Rt, as our radial coordinate 
within the star, Xp = Rp/Ri, » I as our dimensionless pericen- 
tre distance and in = MjM^, as our dimensionless mass. We also 
introduce a fiducial time unit Tq = 2;r(i?p/GMh)''^ and a dimen- 
sionless time T = T/Tq, as well as a fiducial density po = M-„/R\ 
and a dimensionless density p = p/po. The outcome of the disrup- 
tion event depends on the 'penetration factor' y6 = Rp/Rt, that is 
the ratio of the pericenter distance to the tidal radius /?, = q'^^R-^, 
where q = My,/M-^ is the mass ratio between the black hole and 
the star. In order to tidally disrupt the star we require y6 < 1. 
For example, for a mass ratio q = 10*, we have /? = I for 
Rp = lOOR^,. To give an idea of the numbers involved, we note 
that for Rp = lOORg, = lO^Ma, M^, = IMo and R^, = Re, we 
have To » 3.18 10""* yrs w O.Il days, while the unit for the accre- 
tion rate is M^,/To 3.1 lO^Mo/yr. In these units, Eqs. ([!}, ([4} 
and ^ become simply: 



dm 
dr 



-5/3 



3 d6 



(6) 



(7) 



dm f' , , , 

— = 2n I p{x )x dx , 
dx J, 



dm 



dm 
dx ' 



(8) 



(9) 



The above simple set of equations therefore allows us to cal- 
culate the accretion rate onto the black hole as a function of the 
internal stellar structure. In general, we expect the density to show 
a peak at small radii x and therefore a peak at small specific en- 
ergies e. Since material at lower energies contributes to the accre- 
tion at later times, we can already predict what relative changes 
do we expect with respect to the standard r^'' light curve. In par- 
ticular, we expect that if the star is more centrally condensed the 
flare should start with a relatively longer delay (less matter at large 
energies - small return time) and should have a shallower light 
curve (more matter at small energies - large return time). However, 
unless the density is strongly diverging at small radii, we expect 
dm/dx = dm/de to flatten at the lowest energies and therefore the 
light curve to approach a r^'^ profile at late times. 

As an example, we can use the above analytical formulae to 
calculate the specific energy distribution and the accretion rate as 
a function of time predicted for some simple stellar models with 
known density profiles. We have thus considered simple polytropic 
spheres with different indices y = 5/3, 1.4 and 4/3. We have 
first solved numerically the Lane-Emden equation for the three 
cases and have then computed the various relevant quantities us- 
ing Eq. l[6ll-ll9]l above, assuming Xp = 100. The results are shown 
in Fig. [2] The left panel shows the prediction for the energy dis- 
tribution, where the solid line indicates the relatively non compact 
case 7 = 5/3, the short-dashed line indicates y = 1.4 and the long- 
dashed case shows the most compact case y = 4/3. As can be seen, 
the energy distributions do extend up to e ~ 1, but are not flat ex- 
cept at very low energies, the effect becoming more pronounced for 
the more compact cases. The middle panel shows the predicted evo- 
lution of the mass accretion rate m = dm/dr, for the three values of 
7 with the same line styles as the left panel. The red line shows for 
comparison a simple power law with index -5/3. It can be seen that 
indeed the light curves are slightly shallower that r^^^ and approach 
it only at late times. This is even more evident in the right panel, 
where we plot the power law index n = dlnm/dlnr for the three 
cases. If we want to put some numbers on the estimates above, note 
that for our standard numerical values described above a time of 1 
year corresponds to roughly r « 3000. We then see that the power 
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Figure 3. Radial density profiles for the four models considered here. In the left panel we show the density of four solutions of the Lane-Emden equation with 
(from the highest to the lowest central density) y = 1.4, 1.5, 5/3 and 1.8. In the light panel we show the coiTesponding SPH density estimates for the initial 
conditions of our simulations. 



law index after 1 year of the flare is n « -1.5 for y = 5/3, which is 
reasonably close to the expected -5/3. However, such stellar model 
is probably unrealistic for a solar type star, whose structure is rather 
more similar to a y = 4/3 polytrope, in which case, after 1 year of 
the flare the power law index is still n « -0.8. 



3 NUMERICAL SIMULATIONS 

The model described in the previous section is only approximate 
in that it treats the interaction between the star and the black hole 
as instantaneous. In particular, the distribution of specific energy 
of the disrupted stellar material, and consequently the resultant 
lightcurve, has been computed by assuming that the stellar structure 
is essentially unchanged until it reaches pericentre. Still, it high- 
lights some important features of the stellar disruption process: the 
expected energy distribution is in general not flat, and it tends to 
become progressively more peaked towards lower energies as the 
stellar structure model gets more centrally concentrated (that is, as 
the polytropic index y becomes smaller). In order to gain a better 
understanding of the process, we have therefore compared the an- 
alytical expectations with the results of numerical hydrodynamical 
simulations of the process. 



3.1 Numerical setup 

In the case where the encounter is parabolic, as mentioned above, 
the two relevant dynamical parameters are the mass ratio between 
the star and the black hole, q = M],/M-i, and the penetration factor 
p = Rp/Ri- In this work we have considered the case where j8 = 1 
and q = 10'', which imply that the pericentre distance is equal to 
100 times the radius of the star. 

Following the several investigations summarized in the Intro- 
duction, we have also used a non-relativistic SPH code to sim- 



ulate the encoun ter. Our code uses individual particle timesteps 
teate etalJ[T99^ . it evolves the smoothing length by keeping a 
fixed mass within a smoothing sphere (equivalent to roughly 60 
particles) and includes the relevant terms needed to ensure energy 
conservation when the smoothing length is variable (see Price 20051 
for a r ecent review). We also adopt a standard SPH artificial vis- 
cosity jMonaghanl 1 1 992h with viscosity parameters cyspi, = 1 and 
ySsph = 2. 

In order to describe the basic dynamics of the encounter we do 
not require to use an extremely large nu mber of particles in order to 
reach a satisfactory resolution. Indeed, lEvans & Kochane^ ( Il989h 
have shown that their results were numerically converged with a 
number of particles N equal to a few 10*. Even recent calculations 
have o nly used a relatively small num ber of particles, of the o rder 
of 10^ jAval et al. 2000) up to 2 10* jBogdanovic et alj|2004) . In 
this work we have run all our simulations at the two resolution of 
/V = 10* and /V = 10^ and have noticed no appreciable difference 
in the results, thereby confirming the numerical convergence of the 
results. In the following, we only show the higher resolution results. 

We initialize our simulations by placing the SPH particles to 
form the structure of a polytropic star of given index y (we have 
considered the four cases y = 1.4, 1.5, 5/3 and 1.8). This is done by 
initially placing the particles using close sphere packing and then 
differentially stretching their radial position to achieve the desired 
density profile. This method minimizes the statistical noise associ- 
ated with random placing of the particles (we thank Walter Den- 
hen for providing this setup routine). We then relax the structure of 
the star by evolving it in isolation until its internal properties settle 
down. 

We have considered four different values of y = 1.4, 1.5, 5/3 
and 1.8. In this way we encompass the expected range for differ- 
ent kinds of stars, from radiative to convective ones. Indeed, a so- 
lar type star has a density profile close to a y = 4/3 polytrope (it 
is actually best described by y a 1.3). Unfortunately, a y = 4/3 
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Figure 4. (a) Projected density of the of the star at pericentre. The black hole is outside the image, at the origin of the coordinate system. The four panels refer 
to different values of 7 = 1.4 (upper left), 1.5 (upper right), 5/3 (lower left) and 1.8 (lower right), (b) Same as panel (a), but after the encounter, when the star 
is located at roughly two times the pericentre distance. 



polytrope is difficult to simulate, as it has zero binding energy. The 
lowest value of 7 that we use is then 1.4. Red giants and low mass 
stars can be described by a 7 = 5/3 polytrope, while neutron stars 
have a structure which is probably closer to a 7 = 2.5 polytrope. 

We plot in Fig. [3] the initial density profile of our four mod- 
els as predicted from the solution of the Lane-Emden equation (left 
panel) and as realized after the initial conditions have been allowed 
to relax (right). As can be seen, the four models differ in their cen- 
tral concentration, such that the 7 = 1.4 model is the most concen- 
trated and the 7 = 1.8 is the least. It might be worth also to recall 
that models with larger 7 are less compressible than models with 
lower 7. 

Finally, we introduce the black hole as a point mass at the ori- 
gin and we displace the star so as to place its center of mass on the 
required parabolic orbit (since the star is an extended object this ac- 
tually means that the total mechanical energy of the star is slightly 
negative, amounting to roughly -0.005 in our units). The initial dis- 
tance from the black hole is three times the pericentre distance (in 
other simulations not described here, we have also used a larger ini- 
tial distance and found no significant difference). Our code units are 
for length and M^, for mass, which ensure that our results are de- 
scribed in the same dimensionless variables as described in Section 
2. The black hole is modelled as a sink onto which SPH particles 
can be accreted if they come closer to the black hole that a distance 
0.25 in code units. However, in practice, given that our pericentre 
is very large and that we do not follow the evolution of the debris 
long after the interaction, no particles are actually accreted during 
the course of our simulations. 



4 RESULTS 

4.1 The 7 = 5/3 case 

Before comparing the results obtained with various polytropic in- 
dices, we start by describing the results that we have obtained in 
the 7 = 5/3, which is directly comparable to the simulations dis- 
cussed in previous papers. In particular, this simulation is essen- 
tially a higher reso l ution version of the one initially discussed in 
lEvans & Kochanekl ( fT989h . 

Two snapshots of the integrated density profile of the star are 
shown in the lower left panels of Fig.|4la,b), at two different times. 



that is when the star is at pericentre and when it is at roughly 
two times the pericentre distance, after the encounter. The over- 
all structure of the star looks qualitatively similar to the one shown 
in lEvans & Kochanekl ( Il989l) . It is interesting to notice that at peri- 
centre the star is already quite distorted with respect to its initial 
configuration and in particular it has expanded somewhat (recall 
that its initial radius is 1 in code units). This occurs because, in 
isolation, the star is in hydrostatic equilibrium between its pressure 
and its self-gravity. As the star approaches the black hole the tidal 
field effectively acts as to reduce the stellar gravity, making pres- 
sure forces unbalanced and therefore 'inflating' the star. This effect 
is expected to be more significant for small than for large y. This 
reflects the fact that the radius of a polytrope with small 7 is more 
sensitive to the effective gravity. 

A more quantitative comparison can be done by looking at the 
distribution of specific energies of the disrupted star. This is shown 
in Fig.|5]at four different times during the simulations: at f = (up- 
per left panel), at pericentre (upper right), and after the encounter, 
when the star is roughly at four times the pericentre distance (lower 
left) and ten times t he pericentre distance (low er right). For ease 
of comparison with lEvans & Kochanekl ( Il989h . only for this plot 
we have used a logarithmic scale for the distribution. It can be 
seen that initially the distribution is very narrow and centered at 
e = 0, which just reflect the fact that the whole star is initially on a 
parabolic orbit. As the star approaches the black hole, the distribu- 
tion becomes wider and indeed approaches the width predicted by 
the simple analysis of Section 2 (which is equal to unity in the units 
adopted here). The lower left panel, in particular, showing the dis- 
tribution at four times the pericentre, compares almost exactly with 
the distribution shown by Evans & Kochanek ( 1989) (their fig. 3), 
confirming that indeed our simulations replicate accurately their 
results. However, one can see that the density distribution keeps 
evolving until the star is at roughly 10 pericentre distances, where 
it finally settles down in the configuration shown in the lower right 
panel of Fig. [5] We thus see that the distribution is characterized by 
a central peak at lower energies, followed by two 'wings' at larger 
energies. The presence of a central peak is expected based on the 
analytical model described above. The wings, on the other hand, 
refer to the stellar material at the surface of the star, which at peri- 
centre is somewhat distorted from its initially spherical shape (as 
can be seen in Fig.|4| lower left panel) and would obviously show 
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Figure 5. Distribution of specific energies for tlie case y = 5/3. Tiie four 
siiow tile distribution for the initial condition (upper left panel), when the 
star is at pericentre ((upper right), at four times the pericentre distance 
(lower left) and at ten times the pericentre distance (lower right). 



Figure 6. Distribution of specific energies for the case y = 5/3. Solid line: 
average distribution at the end of the simulation. Dashed line: predicted 
distribution based on the analytical model, re-normalized to match the peak. 



some discrepancies with respect to the simple 'spherical' model of 
Section 2. 

Fig. [6] (solid line) shows the distribution of specific energies 
averaged over 10 time units, when the stars has reached ~ 20 peri- 
centre distances and the distribution has settled down. This is com- 
pared with the prediction of the analytical model of Section 2 (cf. 
Fig-El left panel), which is shown with a dashed line. Since the pro- 
files are all normalized to 1, in order to compare the shape of the 
distribution at the peak, we have scaled down the analytical profile 
by a factor x 1.6. We thus see that the analytical profile does ap- 
proximately match the shape of the distribution at the peak, except 
for the presence of the wings, indicating the presence of more ma- 
terial at extreme energies than predicted by the model. Note that, 
obviously, the light curve produced in this case does not show the 
standard r^'^ decline, especially at early times. More details on 
the resultant light curves are given in the next Section, where we 
compare the results obtained with different values of the polytropic 
index. 

4.2 Varying the polytropic index 

We now discuss the effects of varying the polytropic index y on the 
structure of the disrupted star. A first comparison can be obtained 
by looking at Fig. where the various panels show the projected 
density of the star at pericentre (left) and at two times the pericentre 
(right) for the four cases considered (from top left to bottom right: 
y = 1.4, 1.5, 5/3 and 1.8). Several interesting differences can be 
already seen from these images. First of all, note that the overall 
expansion of the star is similar in all cases. However, for larger 
values of y the density structure of the star is much more uniform. 
This is particularly evident in the right panel, which refers to well 
after pericentre passage. In the case where y = 1 .4 the high density 
core is compact, with the density in the 'puffy' tidal tails gently 



declining. In contrast, at the opposite extreme of y = 1.8 the high 
density region is more extended and the edge of the tidal tails is 
more clearly defined, revealing a sharper density cut-off at the edge. 
It is also interesting to note the different degree of internal rotation 
induced in the star by the tidal interaction, with the elongated core 
being more aligned with the line joining the star and the black hole 
(at the origin of the coordinate system) for smaller y than for larger 
ones. 

Let us now look at the distribution of specific energies within 
the star for the four different simulations. This is shown in Fig. [7] 
where the solid lines refer to the simulations, averaged over 10 time 
units when the star has reached a distance of roughly 20 times the 
pericentre. Note that in each simulation, as mentioned above, the 
stars are somewhat inflated once they reach pericentre. In order to 
compare the numerical results with the analytical predictions we 
therefore have to take into account this expansion. We have thus 
simply taken the initial equilibrium density as a function of ra- 
dius within the star and re-scaled the radius by a constant factor 
f , thus effectively applying a homologous expansion to the stellar 
structure. We have then calculated the expected energy distribution 
based on Eqs. l[8j and ^ for this 'inflated' profile. The resulting an- 
alytical predictions are then shown in Fig.|7]with a dashed line. The 
expansion factor to match the numerical data is f = 2.5, 2.1, 1.63 
and 1.6 for the four cases y = 1.4, 1.5, 5/3 and 1.8, respectively. 
It interesting to see that this expansion parameter decreases as we 
increase y, reflecting the reduced response to variations in the grav- 
ity field as 7 gets larger. It can be seen that for y = 1.4 our inflated 
polytropic model describes very accurately the outcome of the sim- 
ulation. However, as y increases the results of the simulations start 
to deviate from the model, in particular in the appearance of wings 
in the tail of the distribution. These wings become progressively 
more prominent as y gets larger. In the previous section we have 
already shown that the core of the distribution for y = 5/3 is well 
described by a non inflated model. Essentially, what is happening 
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Figure 7. Specific energy distribution for the four simulations. Upper left: y = 1.4, upper right: y = 1.5, lower left: y = 5/3, lower right: y = 1.8. The solid 
lines are the result of the simulations, while the dashed lines show the distribution expected from the simple analytical theory outlined in Section 2, where in 
each case the initial density profile of the star has been homologously expanded by a factor f = 2.5, 2.1, 1.63 and 1.6 for the four cases y = 1.4, 1.5, 5/3 and 
1.8, respectively. 



here is that as y increases the expansion of the star becomes pro- 
gressively less homologous, with more material being pushed to 
higher energies (in absolute value). Since the expansion velocity is 
significantly supersonic, the only way to transfer energy within the 
star is through shocks, occurring in the tidal tails. As y increases, 
the density profile of the star becomes shallower, and more material 
undergoes shocks in the outer layers of the star, hence increasing 
the appearance of the wings. We estimate qualitatively the impor- 
tance of shocks in our simulations in the following way. For each 
particle / in our simulations we compute the quantity 

I when (V.u),<0 

qi = < c,j (10) 

I otherwise 



where hi is the particle's smoothing length, (V ■ u), is the local di- 
vergence of flow velocity and c^j is the local sound speed. The 
quantity q is therefore non-zero and negative in regions of conver- 
gent flow and shocks occur where \q\> 1. Fig.[8]shows the structure 
of the disrupted star for the y = 5/3 case at f = 16.75. The top panel 
shows the projected density, while the bottom panel shows a ver- 
tical cross section of q. It can be seen that most of the disrupted 
star is expanding, except for the tip of the tidal tails, where there 
is a strong convergent flow, which has indeed \q\ > 1 and therefore 
undergoes a shock. 



To see how does the effect of shocks changes as the polytropic 
index is varied, we also compute the quantity i5msi,ock> that we de- 
fine as the total mass of particles that have \q\ > 1. Fig. |9] shows 
the time evolution of 5mshock for the four simulations with y = 1 .4 
(solid line), y = 1.5 (short-dashed line), 7 = 5/3 (long-dashed line) 
and 7 = 1.8 (dot-dashed line). This plot shows a few interesting 
features. First, we see that as the index y increases, the amount of 
shocked mass increases as well, confirming our expectation that 
more mass is involved with the shocks in the tidal tails. In par- 
ticular, the two simulations with the largest 7, which are the ones 
displaying the more pronounced 'wings' in the energy distribution, 
are also the two in which more mass undergoes shocks. Second, we 
see that shocks appear to occur in a sequence of peaks. The first 
one, common to all simulations, occurs at ? w 4, which corresponds 
to pericenter passage. For the largest values of 7, we then see a 
number of other peaks, which can be interpreted as the manifes- 
tation of strongly n on-linear stellar pulsation s induced by the tidal 
interaction (see also llvanov & Novikovll2001il . The period of these 
oscillation decreases with increasing 7, consistent with the expec- 
tation that the period of the fundamental mode of stellar pulsations 
should vary as (37 - 4)-"^ (e.g.. Cox 1980). 

As mentioned in Section 2, the fact that for larger 7 the energy 
distribution becomes relatively flatter implies that the distribution 
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Figure 9. Amount of material involved in shocks as a function of time for 
the four simulations with 7 = 1.4 (solid line), y = 1.5 (short-dashed Hne), 
y = 5/3 (long-dashed line) and 7=1.8 (dot-dashed line). As the polytropic 
index grows more mass undergoes shocks, producing progressively more 
pronouced wings in the energy distribution (cf. Fig.0. 



Figure 8. Top: Projected density fot the 7 = 5/3 case at f = 16.75. Bottom: 
veilical cross section of the quantity q, defined in Eq. (To). When \q\ > 1 
the gas undergoes a shock. It can be seen that this occurs at the edge of the 
tidal tails. 



of return times should become steeper and more rapidly approach 
the r^l^ profile expected for an exactly flat distribution. The resul- 
tant accretion rate for the four simulations is shown in Fig.llOKleft 
panel), where the solid line refers to 7 = 1.4, the short-dashed line 
to 7 = 1.5, the dot-dashed line to 7 = 5/3 and the long-dashed line 
to 7 = 1.8. To give an idea of the numbers involved we have ploted 
the results in physical units, assuming = IMq and Rt, = Rq 
(note that if the disrupted star is a giant, the time unit is increased 
by a factor (Rgmnt/Re)^^^, which can be several hundreds, then sug- 
gesting that the rise to peak might be observable, and the decay time 
prior to reaching the asymptotic /"^'^ behaviour can be very long). 
It can be indeed be easily seen that only for the largest values of 
y does the light curve follow the r^'^ profile at early time, while 
for lower 7 the profile gets significantly more shallow. This is seen 
even better in the right panel of Fig. [TO] where we plot the instan- 
taneous power law index n (that is, the logarithmic time derivative 
of the accretion rate) associated with the lightcurve for the cases 
7 = 1.4 (squares) and 7 = 5/3 (triangles) as a function of magni- 
tude drop from the peak (we only plot these two cases for simplic- 
ity: the two other cases follow essentially the same behaviour). The 
dashed line at the bottom indicates n = -5/3. This plot illustrates 
quite clearly that the r''^ regime is only approached at late times, 
after the luminosity has dropped ~ 2 magnitudes from the peak for 
7 = 1.4. In the case 7 = 5/3 the asymptotic regime is approached 
more quickly, after only a luminosity drop of approximately 1 mag- 
nitude. 



5 DISCUSSION AND CONCLUSIONS 

Candidate tidal disruption events of a star by a dormant black hole 
are usually associated with luminous flares in the nucleus of an oth- 
erwise normal galaxy. These c an be detected in X- rays, for exam- 
ple with Chiffldra a nd ROSAT jHalpem et ani2004h or with X MM 
(lEsauei et alj|2008h . or in the optical/UV jOezari et aljlioosh . X- 
ray data generally observe the flare evolving down from the peak 
by a few orders of magnitude, and in the best studied case, NGC 
5905, the decline a ppears to be consistent with a r^^^ fall-off 
jHalpem et alj|2004h . However, in some other cases dOezari et al.l 
12008*) the observations only span a relatively small drop in lumi- 
nosity from the peak. In these cases t he lightcurve appea rs to be 
shallower than and the best fit of lOezari etaP i2008h to their 
optical data indicates a value of n « - 1 . 1 in one case and n « -0.82 
in another. These results are consistent with our prediction that ini- 
tially the lightcurve should be shallow, approaching a r^'^ profile 
only after the luminosity has dropped by 2-3 magnitudes from the 
peak. 

To summarize, in this paper we have revisited the arguments at 
the basis of the expected lightcurve produced by the tidal disruption 
of a star in a parabolic orbit close to a supermassiv e bla ck holes. 
The r ^^^ profile originally proposed bv lReesI ( Il988l) and lPhinnevI 
(Il989h only holds in the case where the energy distribution dm/de 
of the remnant is flat, which we have shown is not the case, in 
general. We have proposed a simple analytical model that relates 
the resultant energy distribution to the density structure of the star. 
This model predicts that more centrally concentrated (solar-type) 
stars should produce flares with a lightcurve shallower than 
approaching it only at late stages. We have tested the model with 
numerical simulations and found that it does reproduce the sim- 
ulated behaviour, with the following two corrections. Firstly, we 
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Figure 10. Left: accretion rate as a function of time for the four simulations. Solid line: y = 1.4 (internal structure close to a solar-type star); short-dashed line: 
y = 1.5; dot-dashed line: y = 5/3 (internal structure appropriate for low-mass, fully convective stars); long-dashed line: y = 1.8. The physical units refer to the 
case where the disrupted star has M-i, = Mq and = Rq (note that if the disrupted star is a giant, the time unit is increased by a factor (Rgiant//?o)^''> which 
can be several hundreds, then suggesting that the rise to peak might be observable, and the decay time prior to reaching the asymptotic r^'^ behaviour can be 
very long). Right: instantaneous power-law index n = dlog M/dlog f for y = 1.4 (squares) and y = 5/3 (triangles) as a function of magnitude drop from the 
peak. The dashed line at the bottom indicates the commonly-invoked power-law index n = -5/3 for the light curve. 



have to account for the inflation of the star from its initial struc- 
ture due to the effective reduction of gravity as it moves in the tidal 
field of the black hole. This is well described by a homologous ex- 
pansion by a factor which becomes smaller as the polytropic index 
becomes larger. Secondly, for large polytropic indices we see the 
appearance of wings in the tails of the energy distribution, indicat- 
ing that some material has been put further away from parabolic 
orbits as a result of shocks in the tidal tails. 

In all cases, we do not obtain a r''^ lightcurve, except at 
late times. Close to the peak of the luminosity, the lightcurve is 
very sensitive to the structure of the star, being shallower for stars 
with polytropic index close to 4/3, expected for solar type stars. 
In this case, the r^'^ profile is reached only after the luminosity 
has dropped by at least two magnitudes. For stars with a relatively 
flat density profile, such as red giants and low mass stars, the r''^ 
profiles is reached earlier. 

In this paper we have only investigated a very simple setup, 
with a given mass ratio between the star and the supermassive black 
hole, and one given set of orbital parameters. It is expected that the 
results would be further dependent on the such additional parame- 
ters, such as the ratio of tidal radius to pericentre distance and the 
eccentricity of the orbit. We plan to consider these effects in subse- 
quent investigations. 

Finally, it should be further emphasised that all these results 
refer essentially to the return time of the disrupted debris, and only 
correspond to an actual luminosity under the further assumption 
that the subsequent accretion is perfectly efficient and oc curs on a 
much shorter timescale, which may not be the case (see iAval et al.l 
I2OOOI) . 
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